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Abstract. From the spectral analysis of a broadband 
(0.1-200 keV) BeppoSAX observation of the Seyfert 1 
galaxy IC 4329A, the main results obtained are: a) the 
amount of reflection, together with the intensity of the 
iron K line, indicate a geometry with a solid angle sub- 
stantially less than 2tt; b) the power law is affected by a 
cut off with e-folding energy about 270 keV, the fourth in- 
dividual object so far where this property has been firmly 
detected; c) two absorption features at about 0.7 and 1 
keV are found, the first corresponding to a blend of O vi 
and O VII, the other to a combination of FeL and NeK 
edges. Compared to an earlier ASCA observation, when 
the source was 30% fainter, the values of both the relative 
amount of reflection and the warm absorber ionization de- 
gree are significantly lower: the comparison is suggestive 
of sizeable delay effects in this object, due to geometrical 
factors in the reflection, and to relaxation to equilibrium 
states in the ionization of the absorber. 

Key words: Galaxies: individual: IC4329A - Galaxies: 
Seyfert - X-rays: galaxies 



1. Introduction 

Starting from EXOSAT and Ginga observations, it has be- 
come progressively clear that the X-ray spectra of Seyfert 
1 galaxies are rather complex (see, e.g., Mushotzky et al. 
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1993, and ref. therein), with the effects of radiative trans- 
fer in the environment of the central engine conrplicat- 
ing the view of the primary emission. The latter is em- 
pirically described as a power law, which, according to 
OSSE/CGRO observations (Gondek et al. 1996, Zdziarski 
et al. 1995), appears to break at high energies in a manner 
that can be described to a first approximation by an expo- 
nential cutoff, although the value of the e-folding energy 
in individual objects remains to be properly constrained. 
At energies below about 1-2 keV the incidence of a fur- 
ther continuum conrponent in excess to the power law, 
the so-called soft X-ray excess, is not yet completely well 
assessed, with ASCA data showing it to be rather less 
common (e.g. Reynolds 1997) than in ROSAT data (e.g. 
Walter & Fink 1993). In fact this issue had a significant 
turn after the discovery with ROSAT and, more accu- 
rately, with ASCA of the rather common incidence of a 
warm absorber: in essence, both types of features are de- 
termined with less ambiguity the better the slope of the 
power law is constrained, a task which ASCA, in addi- 
tion to the higher energy resolution, could achieve more 
directly than ROSAT because of its energy coverage up to 
10 keV. To complicate the matter there is another, seem- 
ingly ubiquitous, continuum component, which affects the 
power law slope determination in the ASCA band (as it 
previously did in the HEAOl and EXOSAT bands) and 
peaks beyond that band, around 30 keV. This component, 
discovered with Ginga, is interpreted as Compton reflec- 
tion of the power law photons off thick matter that can be 
the accretion disk itself and/or the torus envisaged in the 
unified Seyfert 1 and 2 scenario. This component is typi- 
cally accompanied by the iron K fluorescence line, discov- 
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ered in fact before the reflection. The accuracy achieved on 
the strength of the reflection is then relevant to the deter- 
mination of the power law slope and, consequently, of the 
soft excess and/or the warm absorber parameters. Obser- 
vationally the situation is further complicated by the well 
known variability of these objects, such that combining 
data from satellites covering different energy ranges and 
collected at different times leads to further ambiguities. 

The Narrow Field Instruments aboard BeppoSAX of- 
fer for the first time the opportunity to cover simultane- 
ously, and with a reasonably good spectral resolution, the 
band from 0.1 to 200 keV, hence to determine the various 
spectral parameters, especially the broad ones, with better 
confidence than could be done before. This paper reports 
results obtained as part of a broad band spectral survey 
of Seyfert Is, with 2-10 keV fluxes greater than about 1-2 
mCrab. After the observation of NGC 4593 described in 
Guainazzi et al. (1999a) here we present the observation 
of the brightest source in the sample, IC 4329A. In Sect. 2 
the earlier X-ray studies are summarized. In Sects. 3 and 
4, the data reduction and the spectral analysis are respec- 
tively presented. In Sect. 5 these results are commented 
on and compared with previous ones, and the conclusions 
are drawn in Sect. 6. 

2. IC 4329A in X rays 

IC 4329 A is classified as a Seyfert 1.0 galaxy (Whittle 
1992) at z=0.016 (Wilson & Pension 1979). The 2-10 keV 
flux records in the literature show a relatively modest vari- 
ability on timescales of hours-days, and fall within 50% of 
1.0x10"^° erg cm~^ s~^ (corrected for absorption), cor- 
responding to L(2-10 keV) = 10^** erg s~^ (iJo=50 km 
s""'^ Mpc^^). However, from inspection of the 2-10 keV 
light curve from the ASM instrument onboard Rossi XTE 
(quick look results provided by the ASM/RXTE team) it 
appears that this object undergoes occasional outbursts 
by up to a factor of about 10, with timescales of the or- 
der of a few days. Furthermore a factor 3 decline in about 
four days in visible in the light curve from the January 
1993 OSSE/CGRO observation, albeit in the quite differ- 
ent band from 50 to 150 keV (Fabian et al. 1993). 

The spectrum displays neutral absorption at low en- 
ergies, due to a gas column about ten times greater than 
the galactic value N//(,=4.5xl02° cm'^ (Elvis et al. 1989), 
conceivably associated with the disk of the host galaxy, 
which is oriented nearly edge-on (Petre et al. 1984). 

This object is one of the first AGN where a refiection 
component was individually detected with Ginga (July 8- 
10 1989 observation; Piro et al. 1990, Nandra & Pounds 
1994): the amplitude of this component relative to the 
power law was found to be roughly consistent with an 
r2=27r slab of reflecting matter, and the intensity of the 
iron K line simultaneously measured had EW~110 eV 
from a gaussian fit. A time resolved spectral analysis of 
the same data by Fiore et al. (1992) presents tentative 



evidence of a time lag by more than three days in the re- 
sponse of the reflection component to an increase of the 
power law by 20% in three days. Contemporaneous obser- 
vations with ROSAT and OSSE/CGRO in January 1993 
(Madejski et al. 1995) gave for the reflection component 
results, with large errors, consistent with those obtained 
by Ginga, providing in addition a lower limit of about 
100 keV for the e-folding energy of the power law. The 
August 15 1993 ASCA observation has been analyzed in 
several papers (Mushotzky et al. 1995, Cappi et al. 1996, 
Nandra et al. 1997, Reynolds 1997, George et al. 1998). 
The detailed analysis by Cappi et al. (1996), for what con- 
cerns the reflection, provides results surprisingly at vari- 
ance with those obtained with Ginga: despite the flux level 
being only 30% lower, the best fit relative normalization 
of the reflection component is about three times stronger, 
while the gaussian fit of the iron line yields about the same 
EW (89±34 eV at 90% confidence for one interesting pa- 
rameter, i. p.). This result is hard to reconcile with the 
idea that both the line and the reflection arise from the 
same matter, and we shall return to it in Sect. 5. 

At low energies Madejski et al. (1995) discovered with 
ROSAT the existence either of the ionized edge imprint- 
ing by an additional warm absorber, or of a soft excess. 
The ambiguity was convincingly resolved by Cappi et al. 
(1996), who found in the ASCA data evidence of two 
rather strong edges, one consistent with O vii, the other 
with O VIII (this aspect is investigated also in Reynolds 
1997 and George et al. 1998, see Sect. 5). 

3. Observations, data reduction and temporal 
analysis 

A general description of the Italian-Dutch satellite Bep- 
poSAX can be found in Boella et al. (1997a). The ob- 
servations were made with the four, coaligned Narrow 
Field Instruments: the two imaging instruments cover- 
ing, respectively, the band 0.1-10 keV (LEGS, Low En- 
ergy Concentrator Spectrometer, Parmar et al. 1997) and 
the band 1.8-10.5 keV (MEGS, Medium Energy Concen- 
trator Spectrometer, Boella et al. 1997b); the two coUi- 
mated instruments, covering respectively the band 4-120 
keV (HPGSPC, High Pressure Gas Scintiflation Propor- 
tional Counter, Manzo et al. 1997) and the band 13- 
200 keV (PDS, Phoswich Detector System, Frontera et 
al. 1997), with their collimators operating in the rocking 
mode (time-on = time-off = 96 s) to monitor the back- 
ground throughout the observation. The data from the 
HPGSPC, an instrument tuned for spectroscopy of very 
bright sources, provided in our case barely significant con- 
straints and will therefore not be considered in this paper. 
Due to rather subtle problems that remain in the exploita- 
tion for spectral analysis of LEGS data above 4 keV, we 
shall use only data in the range 0.1-4.0 keV. 

In Table 1 the start date of the observation, the net 
exposure times and the net average count rates are given. 
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The exposure time in the LEGS is much shorter than in 
the MEGS because the former instrument is operated only 
during the night-time fraction of each orbit. The reduc- 
tion procedures and screening criteria used to produce 
the hnearized and (between the two MEGS units) equal- 
ized event files are standard and have been described in 
Guainazzi et al. (1999a). In particular, of the two options 
available for the PDS, we used the Rise Time selection 
criterion. 

The spectral counts in the imaging instruments were 
extracted from circular regions of radius 4 arcmin (MEGS) 
and 8 arcmin (LEGS) around the source centroid, and the 
background subtraction was performed using spectra from 
blank sky event files in the same position of the detectors. 
The background count rate used is 2.74xl0~^ c/s in the 
LEGS, 6.95x10-3 c/s (2 units) in the MEGS, and the 
error on the net counts (in Table 1) is dominated by the 
source statistics. Spectra and hght curves from the PDS 
were obtained from direct subtraction of the off- from the 
on-source products, and the error on the net counts (in 
Table 1) is dominated by the background statistics. 

In the field of view of the PDS, 1.3° FWHM, care must 
be taken of sources that could contaminate the target sig- 
nal. IG 4329A belongs to the cluster of galaxies A 3574, 
which is a comparatively very weak X-ray emitter (Pierre 
et al. 1994). We integrated its surface brightness in the 
ROSAT pointed observation described in Madejski et al. 
(1995) to obtain a rough estimate of the flux: in the 2-10 
keV range (with an hypothetically large value of the tem- 
perature, kT=7 keV) it is about 30 times lower than that 
of IG 4329A, hence its contribution in the PDS is totally 
negligible. The same conclusion holds for the point source 
named S3 by Madejski et al. (1995), which lies at about 
12 arcmin from the target, and that we find in the MEGS 
image with a 2-10 keV flux of about 8.4x10"^"^ erg cm~^ 
s~^, more than 100 times fainter than IG 4329A. Another 
source, associated with the elliptical galaxy IG 4329 lo- 
cated only 3 arcmin away from IG 4329A, might contam- 
inate the signal of our target extracted from the LEGS. 
Madejski et al. (1995) attribute to this source a thermal 
plasma spectrum with fcr=0.87keV, and a flux 6.5x10"^^ 
erg cm~^ s~^ in the 0.1-2 keV range: from Fig. 6 of their 
paper we expect the contamination by this source to be 
possibly relevant only below 0.4 keV (but see Sect. 4). 

The detailed study of the observed time variability 
goes beyond the scope of this paper, but some informa- 
tion is presented to justify the spectral analysis of the 
integrated spectra. Fig. 1 illustrates light and hardness 
ratio (HR) curves. The object underwent important in- 
tensity variations, but the x^ test applied to the strings of 
HR values (Table 2) shows that the evidence of deviations 
from a constant value is marginal overall: we conclude that 
spectral variations are unlikely to introduce, in the anal- 
ysis on the integrated counts, an important bias in the 
estimate of the parameters. We shall, though, take care 
of the inhomogeneity in the time coverage between LEGS 



and MEGS, by letting the relative normalization of the 
two instruments as a free parameter. The normalization 
of the PDS to the MEGS wiU instead be held fixed. 



o.g 

0.8 
0.7 
0.6 
1.4 
1.2 
1 
0.8 



If 



o.aa 

0.8 
0.75 

0.7 
0.65 

1.5 



0.8 
0.6 
0.4 



— 1 1 1 1 1 1 1 1 1 1 1 r- 



■++-ui+, 



ECS 0.1-4 keV T "f 






LECS (1.8-4/0.1-1.8) 



++^ 



+-+- 



.++ 



MECS 1.8-10.5 keV 



+ 



++, 



H 1 1 H 



H 1 1 1 1 1 1 1 H 



+ 



t++ V++++++++++ ++ v++^ ++ 



MECS (4-10.5/1.8-4) 



H H 



^ ++++++^++^++^++-^v++^+++ 



PDS 13-100 keV 



PDS/MECS (13-1007l.8-10.5) 



5x10* 10« 

T(s) 



1.5xl0» 



Fig. 1. Light and HR curves of the IC 4329A observation, 
in time bins of 6000 s. (a) Light curve in the 0.1-4.0 keV 
LECS band and (b) HR(1. 8-4.0/0.1-1. 8). (c) Light curve in 
the 1.8-10.5 keV MECS band and (d) HR(4.0-10. 5/1. 8-4.0). 
(d) Light curve in the 13-100 keV PDS band and (f) 
HR(13-100/1. 8-10.5). The LECS and MECS light curves are 
not background subtracted. 



4. Spectral analysis 

The analysis of the integral spectral counts was per- 
formed using the software package XSPEG (version 10). 
In the first step a "Baseline Model Spectrum" (BMS) was 
adopted, which is composed of: a Power Law with an ex- 
ponential cut off (PL), described by three parameters, the 
normalization A, the photon index F and the e-folding 
energy E/ (AxE^^" x exp{—E/Ef))] a Refiection Gom- 
ponent (RG) with two parameters, r ~ r2/27r, the solid 
angle fraction of a neutral, plane parallel slab illuminated 
by the PL photons, and i, its inclination angle to the line 
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Table 1. Observation epoch and mean count rates. The count rate for the MEGS refers to two units 



Start date 


LEGS 


LEGS 


MEGS 


MEGS 


PDS 


PDS 




tcxp 

s 


GR (0.1-10 kcV) 

cts S^"*" 


texp 

S 


CR (1.8-10.5 keV) 
cts s^^ 


tcxp 

S 


GR (13-200 keV) 
cts s^^ 


1998-Jan-02 (9h Im 16s UT) 


25097 


0.696±0.005 


81826 


1.551±0.004 


75214 


0.983±0.015 



Table 2. Hardness ratios: mean values and x^ for a constant value 



HR(cn^)" Xl/A.o.i. 



1.04 



1.27/28 



HR(^^)' Xi/A.o.i. 



0.74 



0.83/28 



HR(ri 



x'/d.o.f. 



0.61 



1.31/28 



"based on LEGS data only based on MEGS data only 



of sight (module pexrav for PL and RC together); a uni- 
form neutral column of gas in photoelectric absorption, 
Nj:/, in addition to the adopted value of the galactic ^Hg 
(module wabs; in both the slab and the column the ele- 
ment abundances used arc the cosmic values in Anders & 
Grevesse 1993); a gaussian iron K line, with three param- 
eters, 'Ejk, cTk aud the intensity I^, also given as an EW. 
Following the procedure adopted by Cappi et al. (1996), 
we included in the model also two absorption edges, with 
four parameters, their energies, Ei and E2, and maximum 
optical depths, ti and T2- 

The energy bins chosen represent about one third of 
the instrumental resolution, which is a function of the 
energy. The normalization C of the PDS relative to the 
MEGS was adopted equal to 0.8 (Fiore et al. 1999). The 
"statistical" errors have been obtained holding G fixed at 
this value, and correspond to the 90% confidence inter- 
val for two i. p., or Ax^=4.61. The "systematic" errors 
linked to the current ±10% uncertainty on G will also be 
investigated. 

The energies of the iron line and of the absorption 
edges are reported as in the rest frame of the host galaxy. 



Table 3. Baseline Model Spectrum fit with cos i=l 



F(2-10 keV)" 


13.17±0.04 


A" 


4.16 


r 


, 0^+0.03 (+0.03) 
J--O"_o.03 (-0.02) 
O-I+0.26 (+0.22) 
"'•"'-'--0.50 (-0.08) 


Nh (10^1 cm-2)<= 


E/ (keV) 


r,y^ + 167 ( + 70) 
^"-•-SO (-40) 


r 


„ t-r+0.15 (+0.27) 
'-'•"'°-0.13 (-0.20) 


Efe (keV)'* 


D-48_o.l6 


ak (koV) 


OMtlil 


Ifc (10^^ cm-2 s-"-) 


.. r„+0.73 (+0.12) 
-^-•""-0.59 (-0.21) 


EWfe (eV) 


if,Q+50 (+8) 


El (keV)"^ 


o.rstlTo 


Tl 


0.52+°;« 


E2 (keV)'' 


i-03t°:i^ 


7-2 


0.19lf,« 


xVd.o.f. 


159.9/141 



"As observed in 10 ^^ erg cm '^ s ^ 

''In 10"^ cm'2 g-i j^gy-i 

'^In addition to NHa=4.55xlO^" cm"^ 

''in the source frame 



4.1. The BMS fit 



The BMS fit was performed on the spectral counts from 
the three instruments together. The results are given in 
Table 3 and illustrated in Fig. 2. According to the x^ 
statistics, the BMS lies above an acceptability limit of 
10%. Moreover, if the contribution to the x^ from the 
LEGS points in the interval common to the MEGS, 1.8-4 
keV, is subtracted, it improves dramatically to 115.9/118. 
We see no evidence in the residuals of a significant contri- 
bution by IG 4329 below 0.5 keV. 



Due to the weak dependence of the RG shape on cos i, 
the inclination angle is very poorly constrained, and con- 
sequently also the value of r, as shown in the confidence 
contour plots in Fig. 3: the best fit value of the angle is 
practically 0°, while the 90% upper limit read from the 
graph is about 70°. Lacking an objective independent es- 
timate of this angle (but see Sect. 4.2), we adopt as a 
reference the face-on configuration for the best fit BMS 
results, as given in Table 3. 
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Fig. 2. LEGS, MEGS and PDS spectra with the best fit 
Baseline Model Spectrum (upper panel), and residuals (lower 
panel) . 



Despite the uncertainty on the angle, the evidence of a 
RC is remarkably clear. However, the skewed shape of the 
confidence contours for the couple of parameters r and F 
(Fig. 4) shows that, despite the ample spectral coverage, 
the two parameters remain statistically correlated, and 
that the fractional uncertainty on r is more sensitive to 
the correlation than the one on F. 




1.86 1.8 

Photon Index 




Fig. 3. Gonfidence (67, 90 and 99%) contours of the parame- 
ters r and cos i, when the inclination angle is a free parameter 
in the BMS fit. 



Fig. 4. Gonfidence (67, 90 and 99%) contours of the parame- 
ters r and r, in the BMS fit with cos i set equal to 1. 



The next remarkable result is the evidence of the expo- 
nential cut-off, with rather tight constraints on the energy 
Elf. This is illustrated in Fig. 5, with the confidence con- 
tours of the two parameters F and E/. When the BMS 
is fitted without the cut-off (Fig. 6), the residuals above 
50 keV show a clear curvature, and the x^ is much larger, 
Ax^=40: according to the F-test, this difference implies a 
probability less than 10^'^ that the cut-off is due to ran- 
dom fluctuations in the counts. Notably the two broad 
band parameters F and r turn out practically the same in 
the two fits, thus showing that the estimate of Ej is quite 
well and independently constrained. 

The confidence contours for the energy and width of 
the iron K line in Fig. 7 show that the line energy is con- 
sistent with the neutral value and that the line appears 
resolved (see Sect. 4.2). 

The detection of two edges at about 0.7 and f keV 
confirms the presence of a warm absorber, whose physical 
state will be further analyzed in Sect. 4.3. 

Finally we consider the systematic errors associated 
with the current uncertainty on the value of the PDS to 
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Fig. 5. Confidence (67, 90 and 99%) contours of the parame- 
ters r and E/, in tlie BMS fit witli cos i set equal to 1. 
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Fig. 7. Confidence (67, 90 and 99%) contours of the width 
((jfe) and energy (Efe) of the iron line, in the BMS fit with cos i 
set equal to 1. The line energy is given in the frame of the host 
galaxy. 
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Fig. 6. LECS, MECS and PDS spectra with the best fit BMS 
without the exponential cutoff (upper panel), and ratio of data 
to model (lower panel). 



MECS relative normalization, C, mentioned above. The 
parameters which are significantly affected by this uncer- 
tainty are Nh^ T, Ey, r and the iron line intensity and 
EW. For these parameters, in Table 3, in addition to the 
statistical errors from the BMS fit, we give in parenthe- 
ses the systematic error estimated from the best fit values 
obtained when C is set at 0.72 and 0.88. Notably, while 
for Nh, E/, Ifc and EW the systematic errors are compar- 
atively small, this is not the case for r, where these errors 
are larger than the statistical errors, and for F, where they 
are comparable. 

J^.2. Further modelling of the iron K line 

The line width is at least three times larger than in the 
optical emission lines (Wandel et al. 1999), and is sugges- 
tive of an origin from the accretion disk. In this case it is 
more natural to replace the gaussian with the relativistic 
profile, as described in the module DISKLINE in XSPEC, 
where a Schwarzschild black hole is adopted. In doing that 
we fixed the line emission energy at 6.4 keV, the innermost 
radius at 6 times Kg=GM/c^, i.e. the innermost stable 
orbit, the exponent of the radial dependence of the emis- 
sivity at /3=-2 (the use of somewhat different values, such 
as /3=-2.5, the average inferred by Nandra et al. 1997, on 
a large sample of objects, leads to very similar results). 
The free parameters are the intensity Ife, the outer radius, 
RoMt in units of Rg, of the disk region contributing to the 
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line and the inclination angle; under the assumption that 
most of the RC comes from the same matter, we required 
the angle to be the same for the two spectral components, 
hence this modification of the model leaves the number 
of free parameters unchanged with respect to the BMS fit 
with the angle fixed at 0° given in Table 3. 

The fit with this model gives a x^ identical to the fit 
with the BMS, and the inclination angle turns out to be 
491^^^ (notably smaller than that of the galaxian disk). 
The line parameters are Rout >158, Ifc=(1.63J;Q;5g)xlO~'* 
cm~^ s^^, corresponding to EWfe^llSJ^Ig eV. Compared 
to the BMS fit results, the values of the other parameters 
in common are practically the same, except for the best 
fit value of r, which, as a consequence of the inclination 
angle being much larger than 0° , turns out to be somewhat 
greater, r=0.72+o;^^. 

One further fit comprises the addition of a narrow line 
at 6.4 keV to the previous model, with the intensity of 
this line being then the only additional parameter. This 
addition leads to a marginally significant decrease in chi 
square, Ax^=-3.7 (93% significance, F~test), with the de- 
tection of a weak narrow line of KW kn— 4:^-31 eV. The pa- 



103;^^ 



rameters of the diskline become Kout <43, EWfe 
eV, the inclination angle i~45°_Q , the reflection nor- 
malization r=0.72lQ^g, the others remain practically un- 
changed (the upper 90% confidence bound on r is much 
smaller than in the previous fit because of the substantial 
reduction of the corresponding bound on the angle) . 

If the gas with the observed Nh completely surrounded 
the source, the EW of the associated narrow fluorescence 
line would be lower than observed by a factor about 10 and 
in order to match the observed value the iron abundance 
should be higher than normal by about the same factor. 
It is more plausible that the narrow line is produced in 
thicker neutral matter at the outskirts or outside the ac- 
cretion disk, which should contribute also to the reflection 
continuum. If we assume that the iron abundance is the 
same in the disk and in this matter, then this contribution 
can be inferred from the ratio of the narrow to the total 
intensity, to be about 1/3. 



luminosity L in the interval 5 eV to 20 keV. We imposed 
the spectral slope of this source to be equal to that of 
the PL. The results of the fit are given in Table 4. We 
note the improvement in x^ with respect to the BMS fit, 
Ax^=-2.2 with two less parameters, while the parameters 
in common, in particular I^ and r, have remained practi- 
cally unchanged. This means that, although a single zone 
ionization model may often be an oversimplification (see 
e.g. Reynolds 1997 and ref. therein), in this observation, 
perhaps due to the limited statistics in the LEGS data as 
compared to ASCA data, it seems to apply rather well. It 
must be recognized, though, that ^ turns out to be poorly 
constrained, likely because of the reciprocal interference of 
the cold and warm gas columns, which are of comparable 
thickness. In Fig. 8 the best fit model, deprived for clarity 
of the contribution by the cold absorbers Nh and ^Hg, 
illustrates the presence of two absorption features corre- 
sponding to the two edges detected in the BMS fit, namely 
a blend of O vi (0.67 keV) and O vii (0.74 keV) edges, 
and a shallow trough around 1 keV, due to contributions 
by Fe L and Ne K edges. 



Table 4. Fit with a warm absorber (cos i=l) 
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4-. 3. The warm absorber 

In order to quantify, in the frame of a single zone, equi- 
librium ionization model, the column N\y and the ion- 
ization parameter ^ of the warm absorber, we fitted the 
data with a model where the four parameters describing 
the two edges in the BMS are replaced by the two pa- 
rameters Nw and £,, and the XSPEC module used is AB- 
SORI. In this module the opacity of the gas is based on 
the ionization distribution of the relevant atomic species 
(same abundances as in the module WABS), in a slab of 
Thomson-thin gas, as a function of the ionization param- 
eter (^—L/nB? (erg cm s^^), where n is the number density 
of the gas and R its distance from the ionizing source with 



5. Comments and comparison with previous re- 
sults 

We now comment on our results, also in comparison with 
earlier results as described in papers mentioned in Sect. 2. 
If we compare the BMS fit with i = 0° (Table 3) with 
the DISKLINE fits (i=45°-50°), we note that the EW of 
the iron line turns out to be practically the same, about 
110 eV. Adopting the "normal" abundances in Anders & 
Grevesse (1993), for a.n Vt = 2tt geometry the expected 
EW is 192 eV with i = 0° (Matt et al. 1997), or about 165 
eV with i = 45°-50° (Matt et al. 1991). The difference in 
both cases with respect to the measured value is significant 
and indicates that fi is smaller than 27r. This indication 
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Fig. 8. The best fit model when the two edges in the BMS are 
replaced by a column of "warm" gas. In this figure the cold 
absorber contribution is ignored for clarity. 



is well supported by the estimate of r: this is significantly 
smaller than one in the BMS fit; in the DISKLINE fit it is 
only marginally smaller than one, but the contribution to 
the RC by material giving rise to the narrow line could ac- 
count for about one third of it, as noted in Sect. 4.2, thus 
bringing the values of r and EW from the disk in closer 
agreement. An i7=27r geometry could be retained by as- 
suming either anisotropy in the primary emission pattern, 
or time delay in the response of fluorescence and reflection 
to changes in the continuum, as tentatively found by Fiore 
et al. (1992). 

The average flux in our observation is nearly identi- 
cal to (30% higher than) the Ginga (ASCA) observations. 
Comparison with our BMS fit of the similar one performed 
by Fiore et al. (1992) on the Ginga observation shows a 
substantial agreement. For a comparison with the ASCA 
results discussed in Cappi et al. (1996), it is appropri- 
ate to report first the outcome of a reanalysis of their 
best fit model, which is equivalent to our BMS with 2=0°, 
performed with the same tools used by us and the 90% 
confidence errors for two i.p.: r=1.98to;!j9, r=2.58ti;n, 
Hue EW=87l42^ eV. We note that both Gamma and fine 
EW are within the errors consistent with our results, but 
that there seems to exist an irreducible discrepancy in the 
value of r. It is then to be remarked that, within the ample 
margins of the error on the iron line EW, the evidence of 
an internal discrepancy with respect to r, noted in Sect. 
2, is by no means compelling. Rather, the comparison be- 



tween the two observations indicates that the intensity of 
the RC can change both in absolute and in relative terms. 
The relative changes could be naturally attributed to a 
delay in the response of the RC to variations in the inten- 
sity of the PL, as suggested also by Cappi et al. (1996) to 
explain the exceptionally large best fit value of r. Their 
suggestion is made the more plausible by the strong, oc- 
casional outbursts of short duration in the ASM/RXTE 
light curve quoted in Sect. 2. 

Concerning the energy of the line, we do not confirm 
the ASCA evidence in Cappi et al. (1996) that it is red- 
shifted by about 100 eV, their gaussian fit value being 
6.30±0.07, but with the errors in our estimate (6.48±0.17) 
we cannot exclude it either. This difference in the gaus- 
sian best fit energy, whatever the cause might be, is the 
most likely cause of the discrepancy in inclination angle, 
from the DISKLINE fit, between us and Cappi et al. (1996), 
who find an upper limit of 25° (same as found also by 
Mushotzky et al. 1995 and by Nandra et al. 1997 with the 
same observation). 

In the ASCA observation Cappi et al. (1996) found 
clear evidence of two edges coinciding with O vii and 
O VIII redshifted by about 25 cV. A similar result, al- 
beit different in the estimate of the optical depths, was 
obtained by Reynolds (1997), who kept the two energies 
fixed at their atomic values. In the BeppoSAX observa- 
tion the ionization degree of the absorber was apparently 
lower, with no significant sign of the O viii edge. This 
difference is further borne out by the estimate of the ion- 
ization parameter. Both Cappi et al. (1996) and Reynolds 
(1997) found £, equal to about 10: we find instead ^ about 
2.5 (although admittedly with a large uncertainty), while 
the best fit values of Nw are comparable. Since the in- 
tensity of the source in our observation was 30% higher 
than in the other, if the strength of the ionizing flux is 
the driving parameter and is proportional to the intensity 
of the power law, the reason of the difference between the 
two epochs must be sought elsewhere. In this respect it 
is remarkable that both Cappi et al. (1996) and Reynolds 
(1997) obtain a worse x^ with the warm absorber model 
than with the two edges model. In addition George et al. 
(1998), which apply to the same observation also more so- 
phisticated models, do not find a single case which passes 
their acceptability limit based on the x^ statistics. This 
is suggestive of a situation which is either more complex 
than the single-zone or far from equilibrium. We propose 
the possibility that at the epoch of the ASCA observa- 
tion the ionization degree was above average and out of 
equilibrium, due to a sharp transition from a high state 
which must have occurred before the start of the obser- 
vation. Non-equilibrium states are discussed in Nicastro 
et al. (1999a), which typically require rather sharp transi- 
tions in the intensity of the source. This possibility is very 
attractive, in that it offers simultaneously an explanation 
for the exceptionally high value of the relative strength of 
the reflection component in the ASCA data. 
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Lastly we stress that we have been able for the first 
time to constrain significantly the energy of the exponen- 
tial cut-off in the PL. The value we find is close to the 
lower end of the ample range estimated by Madejski et al. 
(1995) by combining the simultaneous ROSAT and OSSE 
(Compton GRO) observation with the best matching (in 
flux) fraction of the Ginga observation. That estimate, 
though, is affected by the further uncertainty associated 
with possible variations in the relative intensity of the re- 
flection discussed above, so is the physical modelling of 
the same combination of data by Zdziarski et al. (1994). 
A discussion of the physical implications of our measure- 
ment in terms of self-consistent Comptonization models 
is in progress. 

6. Conclusions 

The main conclusions we draw from our analysis of this 
observation of IC4329A, are the following. 

a) For the flrst time the amplitude of the reflection 
component in this galaxy has been constrained with suf- 
ficient accuracy to draw conclusion from the simultane- 
ously measured, and even better constrained, EW of the 
iron fluorescence line. The values of the two parameters 
indicate that the geometrical factors entering the matter 
illumination by the primary, power law photons imply a 
solid angle substantially less than 27r. A 27r geometry can 
be retained if either anisotropy in the primary emission or 
geometrical lags in the response of the reprocessed pho- 
tons to changes in the primary radiation are deemed likely 
to occur. 

b) The power law is definitely affected by a high en- 
ergy turnover, which, if modelled as an exponential cutoff, 
implies an E/ energy of about 270 keV. At present, this is 
the fourth firm, individual determination of a high energy 
turnover in Seyfert 1 spectra, after NGC 4151 (Jourdain 
et al. 1992, Zdziarski et al. 1996, Piro et al. 1998), MCG- 
6-30-15 (Guainazzi et al. 1999b) and NGC 5548 (Nicastro 
et al. 1999b). 

c) The state of the warm absorber in this observation 
appears different from that found in the ASCA observa- 
tion, in particular the absence of a significant absorption 
at the O VIII K edge implies a lower ionization parameter, 
which contrasts the fact that the intensity of the ionizing 
continuum, as inferred from the direct primary emission, 
is actually 30% higher. We tentatively argue that at the 
epoch of the ASCA observation the ionization state might 
have been higher than in equilibrium, due to a transition 
of the continuum from a high to a low state immediately 
prior to that epoch. Comfortably this same argument of- 
fers an explanation for the strength of the reflection com- 
ponent at the same epoch, which, albeit ill constrained, 
seems to have been irreducibly higher than found at the 
epoch of our observation. In other words, the compari- 
son between the BeppoSAX and the ASCA observations 
is suggestive of delay effects due to geometrical factors in 



the reflection, and to relaxation to an equilibrium state in 
the ionization of the absorber. 
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